Classical limit of non-Hermit ian quantum dynamics 
— a generalised canonical structure 



Eva-Maria Graefe 1,2 , Michael Honing 1 and Hans Jiirgen 
Korsch 1 

1 FB Physik, TU Kaiserslautern, D— 67653 Kaisorslautern, Germany 

2 School of Mathematics, University of Bristol, Bristol BS8 1TW, United 
Kingdom 

E-mail: e.graefe@bristol.ac.uk, hoeningOrhrk.uni-kl.de, 
korsch@physik .uni-kl . de 

Abstract. We investigate the classical limit of non-Hermitian quantum 
dynamics arising from a coherent state approximation, and show that the resulting 
classical phase space dynamics can be described by generalised "canonical" 
equations of motion, for both conservative and dissipative motion. The dynamical 
equations combine a symplectic flow associated with the Hermitian part of the 
Hamiltonian with a metric gradient flow associated with the anti-Hermitian part 
of the Hamiltonian. We derive this structure of the classical limit of quantum 
systems in the case of a Euclidean phase space geometry. As an example we 
show that the classical dynamics of a damped and driven oscillator can be 
linked to a non-Hermitian quantum system, and investigate the quantum classical 
correspondence. Furthermore, we present an example of an angular momentum 
system whose classical phase space is spherical and show that the generalised 
canonical structure persists for this nontrivial phase space geometry. 
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1. Introduction 

While quantum mechanics traditionally demands the Hamiltonian be Hermitian for 
the description of closed systems, non-Hermitian Hamiltonians have proven useful for 
the description of open systems involving decay, scattering, and transport phenomena 
[THE]. Here one uses Hamiltonians whose eigenvalues have negative imaginary parts 
which leads to a monotonic decrease of the overall probability in time. Although 
in most cases non-Hermitian Hamiltonians are introduced heuristically, they can be 
incorporated in a mathematically satisfactory way starting from the system coupled 
to a continuum of states (see [H[3j[7] and references cited therein). Furthermore, 
non-Hermitian VT symmetric Hamiltonians can have a purely real spectrum in 
certain parameter regions and can be used to define a consistent description of closed 
quantum systems [SrllOj. Non-Hermitian quantum dynamics is sometimes regarded 
as a mere perturbation to Hermitian dynamics, adding an overall decay. However, 
non-Hermitian dynamics is in general very different from the familiar unitary time 
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evolution. Much of this difference is due to the non-orthogonality of the non-Hermitian 
eigenfunctions, which can even degenerate at "exceptional points" in parameter 
space. The various intriguing effects associated with non-Hermitian quantum evolution 
attracted considerable attention recently, both from theoretical and from experimental 
side [TTHT9] . The field is still in an early stage of its development and many generic 
features are far from being fully understood. 

Little is known about the quantum-classical correspondence of systems that are 
open in the above sense. Considerable insights have recently been gained from the 
study of toy models described by subunitary or truncated unitary operators associated 
with open quantum maps, where a fractal Weyl law has been found [20H22] . In the 
general case, however, even the classical analogue of non-Hermitian quantum theories 
has hitherto remained an open question. Complex extensions of classical dynamics 
in terms of complex coordinates and momenta have been studied in the contexts of 
both Hermitian and non-Hermitian theories |23If2l)] . But these systems cannot be 
regarded as representing the classical limits of non-Hermitian quantum systems for 
which momentum and coordinate operators have real expectation values. On the 
other hand, in attempts to quantising classical damped motion, a possible link to 
non-Hermitian Hamiltonians has been observed [27lI29] . but studies in this direction 
have not been substantiated, there is no generic structure available, and none of the 
approaches has been able to provide the desired quantum-classical correspondence in a 
rigorous sense. Recently a generalised classical approximation using algebraic coherent 
states [301 - 132] for a non-Hermitian Bose- Hubbard dimer [33,34] was introduced in [35] 
in the context of Bose-Einstein condensates. 

Here we show that there is a generalised canonical structure that can be derived 
as the classical limit of non-Hermitian quantum dynamics. This structure incorporates 
a metric gradient flow generated by the anti-Hermitian part of the Hamiltonian. It is 
closely related to the canonical formulation of classical damped motion suggested in 
the past [36T - I38] . as well as to the gradient flow appearing in thermal and statistical 
physics and control theory [39,40]. It should be emphasised that we do not start from 
a classical damped system and quantise it. We rather begin with the non-Hermitian 
quantum system and perform a coherent state approximation that is often used to 
define the classical limit in the Hermitian case. The generalised canonical structure 
is derived for arbitrary quantum systems having a Euclidean phase space. The Bose- 
Hubbard system studied in [35], however, can be regarded as an angular momentum 
system and thus the classical phase space is given by the Bloch sphere. Here we show 
that the corresponding classical dynamics can nevertheless be formulated according 
to the proposed generalised Hamiltonian structure. This gives some evidence that the 
proposed structure persists for more general phase space geometries. 

The paper is organised as follows: In section[2]we give a brief introduction to non- 
Hermitian quantum dynamics, before we introduce the classical limit using coherent 
states in section [3] We establish the generalised canonical structure for arbitrary one- 
dimensional quantum systems on a flat phase space, using Glauber coherent states. 
This derivation is analogous for higher dimensions. In section[4]we present a case study 
for a damped oscillator. In section [5] we briefly review the classical approximation [35] 
for the angular momentum system studied in [33., 34! and show that the generalised 
canonical structure can be used to describe the classical dynamics on the spherical 
phase space. We end with a summary and short outlook. 
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2. Non-Hermitian quantum dynamics 

For a general quantum system with Hamiltonian H, not necessarily Hermitian, the 
dynamics of a pure state \^f) is governed by the Schrodinger equation 

ih\fy=H\V). (1) 

We can decompose the Hamiltonian into a Hermitian and an anti-Hermitian part via 
T~L = H — if 1 , with H = W and r — T>, where the damping term T is assumed 
to be non-negative. The special choice of the anti-Hermitian part T depends on 
the dynamics to be modelled, and the physical interpretation as a damping will 
become clear later. Under time evolution all states decay towards the subspace 
spanned by the eigenvectors of % with the smallest decay rate (imaginary part of 
the eigenvalue). Assuming a time independent Hamiltonian with a discrete spectrum, 
the time evolution can be expressed as 

\m) = J2 c ^ iEnt/h \^)- ( 2 ) 

71 

Here 'H\ip n ) = £ n \(p n ), £ rL — E n — ir„ with real E n < E n+ i and non-negative T n . The 
expansion coefficients c„ can be obtained from the initial state by projection onto the 
left eigenstates, i.e. the eigenstates of W ■ If there are no real eigenvalues, every initial 
state approaches the most stable initially populated eigenstate in the long time limit: 

\m) -^coe-We-^/'Vo), (3) 

where we have assumed that this is the (non-degenerate) ground state, and the norm 
approaches 

n(t) = (mm)) — > |c | 2 e- 2r ^(^ ) • (4) 
From the non-unitary state evolution (flj one can immediately deduce a 
generalised Heisenberg equation of motion [2] for the diagonal matrix element of an 
operator A (for simplicity we consider only the case dA/dt = 0): 

ih^-mA\>p) = MAii-itfAU) 

dt 

= (1>\[A,H]\1>)-W\M + \'I>), (5) 
where [ , ] + is the anti-commutator. Taking into account the fact that the norm is 
not conserved but decays according to 

h^ = -2(i>\m, (6) 

the dynamical equation of motion for the expectation value of an operator (A) = 
{ip\A\ip} / reads 

ih^- t (A) = ([A,H])-2iA\ r , (7) 

with the covariance 

A ab = ~ (A)(B) . In particular, for A = H we have 

h±- t (H) = -2tf HT . (8) 

This evolution equation simplifies in the special case T = kH with a non-negative real 
constant k to 

h^- t (H) = -2((H 2 )-{H) 2 ). (9) 
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This equation has been suggested by Gisin [41,42., however, with no connection to the 
non-Hermitian Schrodinger equation (JT|). Instead, it was derived from the nonlinear 
quantum evolution equation 

m|0) = fl-|0)-i«f)-f)|^ = ^)-i[|0)<0| ) f]|^> (io) 

with (A) = (<f)\A\<f)) [4D - H5] , One can easily see that there is a direct relation 
between the nonlinear dissipative quantum evolution equation (|10p and the non- 
Hermitian Schrodinger equation (JlJ . A transformation to a renormalised state vector 
\4>) = \tp)/^n using ^ and ©, i.e. hn — —2(T)n, immediately leads to (JTUJ). 
Note that these nonlinear evolution equations can also be expressed in terms of the 
density operator (see, e.g., [4T1I42] ). where they take the form of a double bracket flow 
frequently investigated in control theory and thermal physics |46H48] . Furthermore, 
they can be extended to mixed states evolving towards thermal equilibrium [49 52 . 
Here, however, we intend to investigate the classical limit of the non-Hermitian 
evolution (for a pure state), a topic only briefly addressed in [41 j - 



3. A generalised classical limit for one-dimensional quantum dynamics 

Often the classical counterpart of a quantum system is obtained in a rather sloppy 
manner by simply "taking the hats off the operators" and identifying the resulting 
quantities with the classical ones. For a one-dimensional quantum system describing 
a particle of mass m with momentum p and position q one can conveniently express 
the Hamiltonian in terms of harmonic oscillator ladder operators 

TYtujq + ip + mcuq — ip 

a = =- , a 1 = — ==. (11) 

\2mhu) v2mhoj 

Coherent states \a) are defined by 

a\a) = a\a) , a = j== (12) 

with q = (a|q|a) and p = (a|p|a), and the minimal uncertainty product AqAq = ft/2, 
which can be associated with a phase space point in the limit ft — > 0. 
Applying the identification 

\[ha — > z = — ^ = — and Vha) — > z* = — -J (13) 
\/2moj \/2mtd 

replacing all operators in the Hamiltonian with their associated c-numbers, and then 
identifying the result with a Hamiltonian function on the classical phase space, one 
obtains the classical counterpart in this formulation. The Heisenberg equations of 
motion for an operator F are implicitly replaced with the canonical equation of motion 
in terms of Poisson brackets for the associated function F{q, p) of the classical canonical 
variables: 

P=^[F,H] -> F=j{F,H} ZtZ .={F,H} q , p , (14) 

with Poisson brackets { , }. Here we made use of the complex formulation of classical 
mechanics in terms of z and z* [53 . This classical approximation can be motivated 
in the following way: For an operator function F the relations 

%F] = ^ and [at,F] = -^ (15) 
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hold [53] . Therefore, we can express the Heisenberg equations of motion for the ladder 
operators in the form 

iha = —r and iha = — — . (16) 
oaf da 

We can than think of the identification (1131) as replacing the quantum operators 

by their expectation values in coherent states, for which expectation values can be 

factorised. This immediately yields the canonical equations of motion for z and z*, 

and the corresponding equations of motion for a function of the canonical variables in 

terms of the Poisson brackets. Thus, the classical approximation can be rephrased 

as the assumption that an initially coherent state stays coherent throughout the 

time evolution. This is analogous to the spirit of the so-called frozen Gaussian 

approximation [55., 56 for the dynamics of a quantum system. Note that the 

approximation becomes exact if the Hamiltonian is linear in the group operators a' a, 

a and a) which can be verified by explicit calculation of the action of the resulting 

time evolution operator on an initially coherent state. In a semiclassical context in 

addition one often symmetrises the appearing ladder operators to account for the zero 

point energy. Here, however, we deal with the classical approximation for h — > 

where the zero point energy vanishes identical and thus, we neglect issues of operator 

ordering. Further, while an additional constant changes the energy values, it leaves 

the dynamics invariant. 

Let us now derive the general structure of the classical equations of motion for a 

non-Hermitian Hamiltonian. For this we have to translate the equations of motion for 

the expectation values of the ladder operators for an arbitrary Hamiltonian to their 

classical counterparts using the coherent state approximation. We can express the 

anti-commutators appearing in the equation of motion, making use of the identity |54j 

dF 3F 
[a,F] + = 2Fa+-^, [a^F] + = 2a^F + — (17) 

for an operator function F = F(a,a^). Inserting expressions (fT7f in the generalised 
Heisenberg equations of motion (0 we find 

= (#>-i((2fa+^>-2(a)(f)) 

=-(f>-i((2« t f + I)-2(at)(f)). 

Applying the coherent state approximation hence yields the desired classical equations 
of motion governed by the complex valued Hamiltonian function % = if— iT, where H 
and r are the Weyl symbols of the operators (i.e. the expectation values in coherent 
states) : 

. . _ ag . ar _ou dh .dr _ m* 

dz* dz* dz* ' dz dz dz 

The evolution equation ([6} for the norm is then approximated by 

Hh = -2(r + r )n, (20) 

with constant Tq (note that F = T(z,z*) depends implicitly on time). This evolution 
equation for the norm in fact goes beyond a purely classical description because it 
depends on h. Furthermore, as will be become clear from the examples below, T tends 
to zero in the long time limit and the constant Tq accounts for the finite asymptotic 
decay rate (HJ. 
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Translated to q and p the equations of motion (|19l) can be expressed as 
. OH I dT , . dH dT 

q = — — and p= — muj—. (21) 

op moj dq dq op 

For the time evolution of a dynamical variable A(q,p) with dA/dt — we therefore 

find 

• , , 8A. OA. 
A(q,p) = —q + —p 

dAdH dAdH 1 8A3T mu dAdT , 2 2) 
dq dp dp dq muj dq dq dp dp 



{A H} q , p 



1 dA dT dA dT 
mui dq dq dp dp 



Equation (|2ip can be represented in a generalised canonical structure. To see this 
we introduce phase space variables of equal dimension via the canonical transformation 

q — >• q/y/muj, p — >• y/mjjjp. (23) 

Then equation (|21[) can be written in the form 



= O^VH - G _1 VI\ (24) 



q 
v 

where V is the phase space gradient operator, Q is the symplectic unit matrix and 
G denotes the phase space metric. Since the considered phase space has Euclidean 
geometry, we thus have: 

jHj -o 1 ) ■ G -(l ?) <»> 

(see appendix |Appcndix A| and [ST] for further details) . 

The dynamical equation (1241) is a combination of a canonical symplectic flow 
generated by the real part H of the Hamiltonian function and a canonical metric 
gradient flow generated by the imaginary part T of the Hamiltonian function. The 
symplectic part evidently gives rise to the familiar Hamiltonian dynamics of classical 
mechanics. The appearance of the gradient flow may at first seem surprising in the 
present context. However, if we recall the fact that the gradient vector with a negative 
sign points in the direction of the steepest descent of the function T we see that this 
part of the dynamics drives the system towards the minimum of T and thus can 
naturally be associated with a damping (see, e.g., [36H40] ). We shall show in section 
[3] that this structure persists for a spherical phase space arising as the classical limit 
of the quantum angular momentum system studied in [351157] . 

Finally it should be pointed out that the frequency w appearing in the classical 
equations of motion, i.e. the frequency of the harmonic oscillator introduced to define 
the coherent states, is still a parameter and can be adjusted, for example to improve 
the quality of the classical approximation. 

In summary we identified a generalised classical canonical structure (|24|) 
associated with non-Hermitian dynamics for a flat phase space. We derived this 
structure from a coherent state approximation and it should be noted, that it does 
not depend on the specific choice of H and T. In the case where both are linear in the 
generators of the Heisenberg-Weyl algebra, a, a) and aa\ and for a coherent initial 
state this approximation becomes exact and the quantum dynamics is fully captured 
by the classical equations of motion. To get further insight into the quantum classical 
correspondence for non-Hermitian systems, we shall present some example studies in 
the following section. 
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4. The damped oscillator 

In this section we will study the dynamics of a damped oscillator starting with a purely 
harmonic case. Here the classical approximation is exact, provided that the state is 
initially coherent. This is a well known correspondence identity for Hermitian systems 
which also extends to the non-Hermitian case. In the subsequent section we present 
some results of a case study of an anharmonic oscillator, where the classical dynamics 
is approximate. 



4-1. The damped harmonic oscillator 

Let us first consider the purely harmonic case 

* = + (26) 

and assume that the non-Hermitian part of the Hamiltonian is 

f = kH a (27) 

with a positive real number k. Thus, the metric gradient part of the resulting dynamics 
(|24l) will drive the system towards the minimum of the energy. The classical equations 
of motion (|2H are 



q = p/m — jq , p = —muj 2 q — jp, (28) 

with 7 = ku. Eliminating the momentum p we obtain 

q + 2 iq +(uj 2 +7 2 )< Z = 0, (29) 

which is the familiar classical equation of motion for a damped harmonic oscillator, 
however, with a frequency ujq = \J uj 2 + 7 2 , which implies that the damping is always 
subcritical. The evolution equation for the damped harmonic oscillator (|29|) has been 
presented by Gisin in [4Tj . 

Here the Hamiltonian "H is linear in the generators a, and a^a + 1/2 of the 
harmonic oscillator algebra. Hence an initially coherent state remains coherent and 
the classical evolution equation for q and p agree exactly with the quantum evolution 
of the expectation values (q) and (p) . 

Note, however, that this is no longer the case if one uses (i) coherent states in 
the classicalisation belonging to a harmonic oscillator with frequency uj' ^ uj] or (ii) a 
different damping operator T, for example, T = kp 2 /2m (physically that means that 
the anti- Hermitian part seeks to minimise the kinetic energy only). In the first case 
one obtains the classical evolution equation 

q + j(— + — )q+(uJ 2 +-f 2 )q = 0, (30) 

and in the second case one finds 

q + 2jq + uj 2 q = 0. (31) 

Both results are now approximations of the true quantum dynamics because the 
Hamiltonian is no longer linear in the generators of the algebra. The same classical 
equation of motion also appears in the dynamics of a Caldirola-Kanai Hamiltonian 
which describes the classical dynamics of a harmonic oscillator whose mass increases 
exponentially in time. This Hamiltonian can immediately be quantised (see the recent 
review [58]) and yields dynamics different from the one investigated here. 
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In what follows, we choose the proper coherent states and T as given in (I27[) . In 
this case, the Hamiltonian can be simply written as 

U = h{uj-\ 1 ){a ] a + 1/2), (32) 

a harmonic oscillator with a complex frequency lo = co—i'y, and many results derived for 
the Hermitian harmonic oscillator are also valid here. For instance the time evolution 
of an initially coherent state |ckq) is given by (see, e.g., [51] ) 

\^(t)) = | aQ ) = e i < 2t/2-|ao| 2 (l-e- 2 - t )/2| a ^ (33) 

with 

D t = e -u5t - 1 , (34) 

whereat = (l+-D t )ao = e _I "*ao satisfies the (classical!) equation of motion id t = uiol. 
The time dependence of the norm goes exponentially to zero for long times, 

n t = (ip(t)\^(t)) = e -7^l«o| 2 (i-o- 2 ^) ^ e -nr* ^ o ; (35) 

whereas for short times it decays linearly according to n t ~ 1 — 7(1 + |ao| 2 )i. Note that 
the time derivative of n t agrees exactly with the classical equation (|20|) for To = kfvy/2, 
which justifies the ad hoc insertion of this correction term. 

4-2. The forced and damped harmonic oscillator 

We now consider one of the most celebrated dynamical systems in textbooks, the 
(classical) damped harmonic oscillator with harmonic driving 

q + 27 q + ujI q = Fo cos Qt . (36) 
Here all trajectories approach the limit cycle 

q(t) =Qcos(Q.t- 5), (37) 
where the amplitude and phase shift are 

arctan — * • (3°) 



^ [loI - n 2 ) 2 + 4 7 2 o 2 ' uj 2 - n 

The corresponding quantum system has been studied only rarely (see, e.g., [441150] ). 
Here we consider the Hamiltonian 

U = f£j(a^ a +1/2) + hfta + hf*^ , (39) 

again with a complex frequency ui and a time dependent driving f t . The closed form 
solution known for real frequency (see, e.g., |54[ § 3.11]) can be extended to the non- 
Hermitian case with complex frequency ui. For convenience we provide the solution 
here for an initially coherent state |ao) : 

\lp{t)) = e -^t/2+A t +B tao + (Ct+D t a )a^ a ^ _ e ~iut/2+A t +B t a a - \ a 1 2 /2+\a , | 2 /2 | Q ( ^ ^qJ 

with at = Ct + (1 + D t )ao- The parameters satisfy the differential equations 

lD = 5(£> + 1), LB = / t (£> + l), iC = uC + f;, iA = f t C (41) 
with initial conditions Aq — Bq = Cq = Dq = 0. The solutions are 

B t = -i f&t'e- m 'f t , , C t = -i [ dt'e^'-VfZ , A t = -i Ut' f t ,C v (42) 
J o Jo Jo 
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and D t is again given by (|34l) . Note that a t , parametrising the coherent state, is a 
solution of the differential equations 

id = 2a + / t *, id* = -2*a* - f u (43) 

which are indeed the classical equation of motion 

•• dn ••* dn * 

do* a a 

for the Hamiltonian % — uia*a + f t a + / t *a*, where a and a* act as canonical 
coordinate and momentum, respectively. For the special case of harmonic driving, 
ft = fo cos fit, one obtains 

fa 



B t = 



u 2 - Q 2 



C t = - 



{(£ + fije-^-")* + (2 - fi)e- i(< " +n)t - 22}, (45) 
- { (2 - fi)e im + (2 + f7)e" i0t - 2we- i5t } , (46) 



/o 



2(2 2 - ft 2 

+ ~2?o2 (( C + ^)e" i(S ^ 0)t + (2 - n)e-^ +n » - 22)}, (47) 

which approach, in the long time limit, to the simpler expressions D t — > — 1, Bt — > 0, 
A t -> A t °° +4 0) , with 



4(2 2 - ft 2 ) 
and 

Ct ^ ~ 2p^) { (5 " + (S + ° )6 " m } = ^ (50) 

with a t — > af, the limit cycle (|37l) . We note that a ( ' t ' > and ^j -* are T-periodic and 
the limiting state for t — > oo can be written in the form 

|Vb(t)> = e- i£oi/a |0o(*)) with \Mt)) = ce-^ 0)+ l Q *l 2 / 2 |a ( t 0) ) = |0 O (* + T)) , (51) 
where c is a constant, and 

-f('-^??)- < 52 > 

The quantum dynamics can be computed numerically, e.g., in the discrete 
operator representation [59] . As an example, figure [1] shows the evolution of a coherent 
state initially located at (qo,po) = (2, 0). The parameters are chosen as uj = 1, m = 1 
(these are kept fixed in all calculations) and 7 = 0.1, h = 1. For the driving force we 
used Q — 1 and fo = 0.1. The time dependence of the expectation value (q) shown as 
a solid black curve agrees, of course, with the analytic formula given above and also 
with the classical oscillation q(t). Asymptotically, the motion approaches the limit 
cycle (03). 

If the initial state is chosen differently from a coherent one, the agreement between 
classical and quantum evolution breaks down, as in the Hermitian case. As an example, 
the dashed blue curve in the figure shows the results for an initial superposition of two 
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Figure 1. Time evolution of the mean position (q) for a damped quantum 
harmonic oscillator (oj = m = ti = 1, 7 = 0.1) which is harmonically driven 
(/o = 0.1, f! = 1). The solid black curve shows the evolution of a coherent state 
initially located at (qo>Po) = (2,0). This curve agrees with the classical motion 
and converges to the limit cycle. The dashed blue curve is obtained for an initial 
superposition of two coherent states located at (2,0) and (—2,0). 



coherent states located at (qo,po) — (2, 0) and (—2, 0) (a so-called cat state). Here (q), 
which is initiallly equal to zero, cannot be described by a single classical trajectory 
but shows an interference structure before it converges to the asymptotic limit state. 
This behaviour is characteristically different from the Hermitian case, where no limit 
cycle exists and the quantum evolution never approaches the classical dynamics. 

Let us finally remark that in the quantum case most states approach 
asymptotically to a unique state, a quantum limit cycle. This is evident, if one 
formulates the dynamics in terms of non-Hermitian Floquet states. Here the time 
evolution can be expressed in the form of a linear combination as in @, where the 
eigenstates are replaced by the Floquet states 



and the eigenvalues by the (complex) quasienergies e n . Then each state approaches 
asymptotically the most stable Floquet state, which can be identified as the quantum 
limit cycle. In the present case of a forced harmonic oscillator the quasienergies are 
known analytically (see, e.g., [591(60] ) and the solution can be directly extended to 
complex frequencies: 



which agrees for n = with the result already given in (f52"j) . The most stable state 
n = is, of course, only reached if the projection of the initial state onto this state is 
different from zero. Otherwise, the state approaches the most stable initially populated 
quasienergy state. In fact, all quasienergy states appear as quantum limit cycles, as 
already remarked by Gisin |44) . 

As an illustration figure [5] shows an example of the quantum limit cycle in phase 
space. Shown in false colors is the Husimi distribution, i.e. the projection on coherent 
states 



\M*)) = e- i£ " t/r V„(i)) with \<p n (t + T)) = \tp n (t)) 



(53) 




71 = 0, 1, ... 



(54) 



P(jp, q; t) — \ia\ij)(t))\ , a — (mojq + ip) / V 2mhuJ 



(55) 
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Figure 2. Husimi distribution of the quantum limit cycle for a damped quantum 
harmonic oscillator (u> = m = H = 1, 7 = 0.1) which is harmonically driven 
(/o = 1, f! = 1), averaged over one driving period. The solid black line shows 
a classical trajectory started at (<7o,Po) = (2,0) for comparison. This curve 
approaches the classical limit cycle that coincides with the ridge of the averaged 
Husimi distribution of the quantum limit cycle. 

which is averaged over one period of the driving in the long time limit. Also shown 
is a classical trajectory started in at po = 0, qo = 2 which converges to the classical 
limit cycle that concides with the ridge of the quantum distribution. 




Figure 3. Comparison of the quantum and classical damped anharmonic 
oscillator with = 0.4 and weak damping, 7 = 0.01. The quantum expectation 
value (q) (solid black curve) is shown as a function of time in comparison with 
the classical evolution q(t) (dashed blue curve). 



Classical limit of non-Hermitian quantum dynamics 12 
4-3. A damped anharmonic oscillator 

As an example of an anharmonic oscillator we consider the Hamiltonian 

1 9 mu! 2 9 B a , 
H = t-P 2 + -^-f + jg 4 , 56 

2m 2 4 

where the non-Hermitian part T of the Hamiltonian is chosen as 

7 ( 1 ^2 



u \2m 2 J 

as for the harmonic oscillator. Then the classical equations of motion (|2"Tj) arc 

q = p/m — jq , p — —mw q — Bq 3 — ■yp. (58) 

In the Hermitian case 7 = one finds the well-known deviations between the quantum 
and classical evolution showing breakdown and revival phenomena. As illustrated in 
figure [3l these phenomena survive in the case of a weak damping, 7 = 0.01, and 
moderate anharmonicity B = 0.4. The oscillation frequency and the overall decay is 
reproduced in the classical approximation. 

Numerical results for the damped case with stronger damping are shown in figure 
[4j where we have used the same method and parameters as described above. The 
damped anharmonic evolution shows differences between the classical and quantum 
evolution and, as expected, the deviations increase with anharmonicity B. This is 
due to the increasing deviation of the quantum state from a coherent one during the 
time evolution. For short times, when the state is still almost coherent, we find a 
very good agreement between classical and quantum evolution. In further numerical 
studies we observed that the qualitative behaviour is not changed, if we also chose the 
anti-Hcrmitian part T anharmonic. A detailed investigation of the behaviour of the 
quantum-classical correspondence in dependence on the special choice of Hermitian 
and anti-Hermitian part as well as their relation is a promising starting point for future 
studies. 



5. An angular momentum system 

Recently a generalised classical approximation in the spirit introduced in the preceding 
sections has been performed for a non-Hcrmitian angular momentum system in the 
context of many-particle mean-field correspondence for ultracold atoms in an open 
double well trap [35,57 . The Hamiltonian investigated there is of the form 

U = 2(L Z + 2vL x + 2cL 2 z - 2i 7 (L z + L) , (59) 

where the Lj are angular momentum operators and L denotes the angular momentum 
quantum number. By Schwinger's harmonic oscillator representation of angular 
momentum, this can be interpreted as a second quantised many-particle model and the 
classical approximation arises as a mean-field approximation. Here we briefly review 
the generalised coherent state approximation for this angular momentum system and 
show that the resulting mean-field dynamics can be described by the generalised 
canonical structure (|24p . however, on a spherical phase space. 

To perform the classical approximation we have to replace expectation values in 
the quantum Heisenberg equations of motion with their values in a coherent state. 
For angular momentum operators the classical approximation will be performed using 
SU(2) coherent states. The classical limit is then realised for large quantum numbers, 
that is, L — > 00, where we set ft = 1 for convenience. 
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Figure 4. Damped harmonic and anharmonic oscillator: Expectation value (q) 
as a function of time for 7 = 0.1 and j3 = (upper left) = 0.2 (upper right), 
/3 = 0.4 (lower left) and /3 = 1.0 (lower right). The quantum expectation values 
(g)(t) (black solid curves) for ft = 1 are compared with the classical evolution q(t) 
(blue dashed curves). 

The dynamics of the angular momentum expectation values follows from the 
generalised Heisenberg equation of motion ([7]): 

&{L X ) = - 2e(L v ) - 2c({L y ,L z ] + ) - 4 7 A|^ 

&(L y )=2e(l x ) + 2c([L x ,L x ] + )-2v(l x )-4r f Al vtm (60) 

±{L z )=2v{L y )-A 1 tf Lziz , 

and the norm of wave function decays according to 

A<*|*> = -4 7 {<4)+£} (61) 

The generalised SU (2) coherent states, often also denoted as atomic coherent states, 
can be constructed by an arbitrary SU (2) rotation of the extremal angular momentum 
state \L) (spin up): 

\6,<t>) =R(6,(t>)\L) =e iHL, S ^-L y cos<t>)^ L y (g 2 ) 

The coherent state approximation — that is, the assumption that an initially coherent 
state stays coherent throughout the time evolution — is exact if the Hamiltonian is 
linear in the generators of the SU (2) algebra, i.e. for c = in our case. This can be 
seen by an explicit calculation of the action of the time evolution operator 

U(i) = exp {-iHt/h) (63) 

on an initially coherent state (1621) . The classical equations of motion are obtained from 
the quantum dynamics of the relevant expectation values by replacing all expectation 
values with their values in coherent states and identify these as the classical quantities. 
The SU(2) expectation values of the Lj, j =x,y,z appear as the components of the 
classical Bloch vector: 



8j = {L 3 )/2L. 



(64) 
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The expectation values of the anti-commutators appearing in (1601) in SU(2) states can 
be shown [3S1[57] to factorise as 

([£i,£ 3 -] + > = 2(1 - l L ){L i ){L j ) + 5 ij L (65) 

Inserting these expressions into (|60|) and taking the classical limit L — > oo (where we 
have to keep 2Lc = g fixed) we obtain the desired non-Hermitian classical evolution 
equations: 

4= -2es y -Ags z s y +47 s z s X7 

s' y = +2es x +Ags z s x -2vs z +4js z s y , (66) 
s z = +2vs y -7 (1 - 4s 2 z ) . 

In the case g = in which the assumption that the many-particle state stays coherent 
in time is exactly fulfilled, these equations exactly coincide with the quantum dynamics 
for arbitrary L for an initially coherent state. The nonlinear non-Hermitian Bloch 
equations are real valued, and conserve s 2 = s 2 + s 2 + s 2 z = 1/4, i.e. the dynamics is 
regular and confined to the Bloch sphere. The decay of the total probability can be 
approximated by 

h = -4 7 L(2s z + l)n (67) 

The classical dynamics is analysed in detail in |35ll5"T] . Here we show that it can be 
described by the generalised canonical structure (|24|) . For this purpose we introduce 
the canonical coordinates p and q on the sphere related to the Bloch coordinates via 
1 



s ^ = ^ V 1 -P 2 cos(2g) 

s y = \ \I^-V 2 sin(2g) 
1 

s z = - 2 P- 

The equations of motion then read: 



q = e + gp-v cos(2g) (69) 

V 1 ~P 



p = - 27(1 - p 2 ) + 2vy/l - p 2 sin(2q). (70) 

This can be expressed as generalised Hamiltonian equations of motion of the form 
where J? is the usual symplectic matrix (|25l) and 

V U 2(l-p 2 ) / 

is the corresponding Kahler metric on the Bloch sphere (see appendix |Appcndix A| 
for some details). The Hamiltonian function T~L = H — IT is given by the SU(2) 
expectation value of the quantum Hamiltonian (|59j): 



H = ep + v\fi -p 2 cos(2q) + |p 2 and T = jp. (72) 

In figure [5] we show an example of the quantum dynamics in comparison with the 
classical evolution on the Bloch sphere for a strong nonlinearity g = 1.5 and a small 
non-Hermiticity 7 = 0.1. Here the classical dynamics has several fixed points, one of 
them being a sink of the dynamics close to the south pole of the sphere and another 
being a source close to the north pole. The initial condition in the present example 
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Figure 5. Quantum (black solid line) and classical (blue dashed line) angular 
momentum evolution on the Bloch sphere. Plotted are the expectation values 
of the quantum angular momentum components in comparison with the classical 
Bloch vector for an initially coherent state centered at the classical initial condition 
close to the north pole of the Bloch sphere up to a time t e nd = 25 for 7 = 0.1 
and g = 1.5 and a total angular momentum of L = 40. 



is close to the source. For short times the behaviour of the quantum and the classical 
dynamics is similar. For longer times, however, while the classical dynamics moves 
towards the sink staying confined on the surface of the Bloch sphere, the quantum 
dynamics can tunnel through the sphere. Finally they both approach the sink of the 
dynamics. Note that the quantum classical correspondence for this system can be 
quite intricate. Depending on the initial conditions and the total angular momentum 
the quantum dynamics can differ considerably from the underlying classical behaviour 
due to interference and tunneling effects. This is analysed in more detail in [57 . 

6. Conclusions and Outlook 

In conclusion we derived a generalised canonical structure (|24p incorporating a 
metric gradient flow that can describe both conservative and dissipative motion as 
a classical limit of non-Hermitian quantum dynamics. The classical limit was defined 
as a coherent state approximation and it was shown that the canonical structure 
arises for Euclidean phase space geometry associated with Glauber coherent states, 
irrespectively of the special choice of the Hamiltonian. Furthermore, we presented an 
example of an angular momentum system where the classical limit is given by the limit 
of large angular momentum. The appropriate coherent states in this context are the 
SU (2) or atomic coherent states which lead to a spherical classical phase space. We 
showed that the resulting classical equations of motion can also be described using 
the generalised canonical structure derived for the Euclidean case. This strongly 
suggests that the proposed structure holds for more general systems with different 
metric structures on the classical phase space. 

The quantum classical correspondence arising from these non-Hermitian systems 
is in many cases equivalent to the well-known behaviour observed in Hermitian 
systems. However, new effects due to the gradient part of the motion are to be 
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expected, the investigation of which is a promising task for future studies. Interesting 
questions here regard the timescales of the quantum classical correspondence (that is 
the non-Hermitian equivalents of the Ehrenfest and Heisenberg times for Hermitian 
systems) and their dependence on the Hermitian and anti-Hermitian part individually 
as well as the relation of Hermitian and anti-Hermitian part. It is also interesting 
how the generated classical flow depends on this relation between the Hermitian and 
anti-Hermitian part of the Hamiltonian and in particular, which simplifications are to 
be expected if they are proportional. 

Furthermore, the discovery of a general structure in the classical limit of non- 
Hcrmitian quantum theories paves the way for the development of new semiclassical 
techniques such as quantisation conditions and trace formulas for these systems. 
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Appendix A. The Metric, the symplectic structure, and a generalised 
canonical evolution on phase spaces 

In this appendix we provide the minimal basic facts concerning the notion of a Kahler 
structure on manifolds that is relevant for the generalised canonical structure of the 
non-Hermitian dynamics. The interested reader is referred to 61,62 for further detail. 

In Riemannian geometry a manifold is equipped with a measure of distance 
determined by the metric tensor g in the following way: The infinitesimal distance 
ds between a point (x , • • • , a;") and (x 1 + dx , ■ ■ ■ ,x n + dx n ) on an n-dimensional 
manifold X n is given by the expression 

ds 2 = g ab dx a dx b . (A.l) 

Here we make use of the Einstein sum convention. The metric tensor thus consists of 
n(n + l)/2 quantities g a b that define the notion of distance on X n . 

Clearly distance is invariant under local coordinate transformation of the form 
(x 1 , • • • , x n ) (it 1 , • • ■ , u n ). In other words, if we write g^ for the metric tensor in 
the transformed (tt , • • • , u n ) coordinates we have 

ds 2 = g]£ du a du b = g ab dx a dx b , (A.2) 

from which it follows that the metric transforms according to 

(u) dx c dx d 

ft? ft?** 

If the manifold in question is a classical phase space then it comes equipped with 
an additional symplectic structure oj. A space accommodating both a metric g and 
a symplectic structure oj is said to possess a Kahler structure if these two quantities 
satisfy the compatibility condition [62], which in matrix form is written 

oj- 1 = (g-'ojg^f, (A.4) 

where the superscript T denotes matrix transpose. It is evident from (|A.4|) that 
the compatibility condition is invariant under the simultaneous scale transformation 
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g — > Kg and u> —¥ kuj for some k / 0. Hence, the Kahler metric and symplectic 
structure are defined only up to a scale factor. On the other hand, if either oj or 
g is given, then this automatically fixes the relevant scale factor according to the 
compatibility condition. 

In classical mechanics one usually associates the symplectic structure on the phase 
space with the canonical equations of motion by demanding them to read 

' I) (A-*) 
and hence fixing the symplectic structure fl as 

fl=(; "o 1 )- (A.6) 

Thus, if the phase space is a Riemannian manifold and therefore equipped with a 
metric g the corresponding Kahler metric G is connected to g up to a factor that is 
determined by the choice (|A.6[) of the symplectic structure. It is straightforward to 
show that for a Euclidean space this coincides with the usual choice of a Euclidean 
metric G = 1. 

While in the conventional context of Hamiltonian dynamics the metric is of no 
consequence, in the present paper we have a complex Hamiltonian function % = H—iT, 
and the metric structure appears as a natural choice to extend the canonical equations 
of motion to the generalised form 



Q^VH - G _1 Vr, (A.7) 



<l 
P 

where Q is given by the usual expression (|A.13[) for Hamiltonian dynamics and G is 
the associated Kahler metric. 

Let us now proceed to calculate the Kahler metric that is compatible with the 
choice (IA.6I) for the symplectic structure, in the case of the Bloch sphere. The 
Riemannian metric g on a two-sphere S 2 C R 3 of radius R, embedded in the three 
dimensional Euclidean phase space via 

x = R sin 9 cos <j> 

y = Rsm9sm<j) (A.8) 
z = RcosO, 

can be found from the transformation rule (|A.3[) (which extends to the case in which 
we have a parametric subspace). Starting from the metric in R 3 given by the identity 
matrix we deduce at once the well-known expression of the metric on the sphere: 

«**-(? *ir,)- (A - 9 > 

Because the sphere is the Bloch sphere the radius is given by R = |. If wc employ 
the usual coordinate system p and q according to 

p = cos9 and q = <j>/2, (A.10) 

that is, 

9 = arccosp and = 2q, (A. 11) 

then we find that the metric of our Bloch sphere is given by 

g M = I f T=F / 2N V (A.12) 
y 4 V 4(1 -p 2 ) ) y ' 
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From this metric we can calculate the symplectic structure oj on the phase space (p, q) 
making use of the Kahler compatibility condition and find: 

«-(° - 2 ). (A.13) 

As a consequence, by comparing (|A.13|) and (IA.6|) we see that the scaling factor n is 
given by k = 4, that is, J? = ^w, and we conclude that the Kahler metric associated 
with the symplectic structure appearing in the canonical equations of motion is given 
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